Conductance through a One-Dimensional Correlated System: 
Relation to Persistent Currents and Role of the Contacts 
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Based on a recent proposal [O.P. Sushkov, Phys. Rev. B 64, 155319 (2001)], we relate the quantum 
conductance through a sample in which electrons are strongly correlated to the persistent current 
of a large ring, composed of the sample and a non-interacting lead. A scaling law in the lead 
length allows to extrapolate to a well-defined value of the conductance, depending only on intrinsic 
properties of the sample and the nature of the contacts between the sample and the lead. For 
strongly disordered samples, the conductance is found to be enhanced by the interaction. 
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Viewing quantum transport as a scattering problem 
[Q, |2j generated a new understanding of the electronic 
conductance. This approach is able to explain a wealth 
of experimental results |^ in mesoscopic systems when 
electron-electron (e-e) correlations are not important. To 
include these correlations is non-trivial and remains one 
of the major challenges in the field (see e.g. While 
none of the proposals to calculate the conductance for 
a correlated system is well suited for numerical calcula- 
tions j|, ^ or free of certain assumptions 0, such an 
issue becomes crucial in present day's research exploring 
electronic transport through nanosystems (carbon nano- 
tubes molecules and point contacts 0), where 
the Coulomb repulsion leads to important correlations. 

Reservoirs and leads are key elements in the scattering 
approach, and possess very clear physical meanings since 
the measurement is made with electrodes which behave 
as electron reservoirs. In a good electrode, the electron 
density Ue is large, the ratio Vs between Coulomb and 
Fermi energy is small, hence the e-e interaction is negli- 
gible. In contrast, rie in a nanosample can be very small, 
yielding a large ratio and important e-e correlations. 

The dimensionless conductance g does not only depend 
on the intrinsic properties of the sample, but also on 
the way it is connected to the electrodes. The quality 
of the contacts is particularly important for correlated 
electrons. For a clean Luttinger liquid attached to non- 
interacting leads through reflectionless contacts, it has 
been found |11 that the interactions do not influence g. 
In the other extreme, if the contacts are tunnel barriers, 
the interactions lead to Coulomb blockade ||l^, thereby 
dominating g. In carbon nanotubes, various transport 
regimes are observed depending on the nature of the con- 
tacts |§. 

As shown by Kohn and Thouless |l^, g is also 
related to the sensitivity of the sample's eigenstates to 
a change of the boundary conditions. This sensitivity 



can be tested by closing a system to a ring and measur- 
ing the persistent current as the response to an enclosed 
magnetic flux (p. At zero temperature, the persistent 
current is given by J = —dE/dcf), where E is the ground 
state energy of the many-body system. Interactions play 
an important role for J, and it is generally accepted that 
they account for the large difference between experiments 
and one-particle calculations |l5| . There have been var- 
ious attempts ^ to link J and g for an interacting 
ring. However, the ring built from the sample itself does 
not contain any reservoirs in which energy relaxation can 
take place. Negative zero-frequency conductivities occur 
|]l7| , unlike in the dissipative case in which we are inter- 
ested here. 

As pointed out in Refs. |5[ |l^, at zero temperature, not 
only for the non-interacting case, but also for correlated 
samples, g is given by |i(i?i?)P, the probability for an 
electron at the Fermi energy Ep to be elastically trans- 
mitted through the sample. Moreover, if one replaces 
the massive electrodes (with negligible e-e correlations) 
used in a real measurement by very long non-interacting 
one-dimensional leads, one can expect that they have a 
similar effect Sushkov recently proposed jlO) that 
\t{Ep)\'^ can be extracted from the persistent current of 
a much larger ring, composed of the sample itself, to- 
gether with a long lead closing the system. This has the 
considerable advantage that a ground state property (J) 
suffices to determine g. However, one needs the J of the 
combined system (sample plus lead), and not the one of 
the system alone as in previous works 0, |l3|, |lj, |l6|| . 

In the following, we adapt the approach of Ref. to 
calculate g for one-dimensional interacting electrons us- 
ing the Density Matrix Renormalization Group (DMRG) 
algorithm ||l^, ^ . We check that a scaling law allows to 
extrapolate to an infinite lead, yielding g as a property 
of the sample and the way it is connected to the lead. 
The nature of the contacts turns out to play a major 
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role. Then, we extend our analysis to disordered sam- 
ples, where we find that, similarly to the case of persis- 
tent currents |2^, repulsive interactions may increase g 
for strong disorder. 

We first present an alternative derivation of Sushkov's 
result 1 10 , pointing out the main assumptions, for the 
non-interacting case. As depicted in the upper inset of 
Fig. 1, we consider a sample (S, hashed region) closed 
to a ring by a non-interacting and disorder- free lead (L), 
and threaded by a flux (j>. The total length L = Lg + Li, 
consists of the sample length Lg and the lead length Ll. 
The one-particle eigenstates of the total system satisfy 

det(/-AfsA/L) = 0, (1) 

with the transfer matrices of the sample and the lead 
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respectively. Here, $ = 27r<l)/(j)Q where 0o is the flux 
quantum. The scattering is characterized by the angle 
8, the phase-shift a, and the angle (p (equal to tt/2 if 
right-left symmetry is respected) . These angles are func- 
tions of fc, the wave- vector in the lead. The transmission 
amplitude is given hy t — e"^" sin 9 sin ip . With (H), the 
quantizing condition (|l|) can then be written as 
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with the relative phase-shift 6a — a — kLg. The per- 
sistent current carried by a one-particle state (with en- 
ergy e) is j{4>) = —{de/dk){d(l)/dk)~^. We work at 
$ = 7r/2 and establish two crucial assumptions: i) 
\d{Sa)/dk\ < L; ii) de/dk ~ h^k/m. The first one 
states that the Wigner time associated with the scat- 
tering region is negligible compared with the time spent 
in the leads. Notice that we work with a relative Wigner 
time T\Y — [m/ h^k)d{5a) / dk, that is, the difference be- 
tween the delay time of the scattering region and that 
of a potential-free region having the same length. The 
second assumption implies that the dispersion relation is 
essentially unaffected by the scattering potential. 

The persistent current of N non-interacting spinless 
fermions (for simplicity we take N even) is given by pl|] 



J{^^^/2)^^kF\t{kF)\- 
mL 
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Denoting by J° the persistent current of a clean ring of 
length L, the conductance may be expressed as llG] 



= lim 



J(7r/2) 
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FIG. 1: Scaling of InD with the total length L of the system, 
showing a linear increase of In D with 1 / L for even Ls = 12 
{U = 1 (•), U = 2 (■), (7 = 3 (A), [7 = 4 (♦)), and a de- 
crease for odd Ls = 13 ((/ = 1 (o) and (7 = 2 (□)). Lower 
inset: scaling for unpolarized electrons within the Hubbard 
model (Ls = 2, (7 = 1). Upper inset: sketch of a ring consist- 
ing of the sample (hashed region) and a non-interacting lead 
threaded by a flux (j). 



With interaction |22 , assumptions i) and ii) always 
hold in the large L limit. Moreover, the use of Eq. ^ 
implies that the one-particle states of the correlated sys- 
tem can still be indexed by the lead wave- vectors k. That 
is, adding an infinite non-interacting lead to a finite non 
Fermi- liquid sample restores the Fermi-liquid behavior. 
This assumption, which has been used in the pertur- 
bative calculation of transport through Hubbard chains 
connected to reservoirs p3| , requires that the interactions 
are completely switched off in the one dimensional lead. 
Otherwise the Luttinger liquid behavior p4) sets in, and 
one cannot obtain a result which is independent of the 
length of the auxiliary lead. In this, our approach differs 
from Sushkov's, where the interactions in the lead are 
kept (within the Hartree-Fock approximation). 

Eq. (||) allows to calculate g from ground-state ener- 
gies. We do this for spinless fermions (polarized elec- 
trons) in a ring described by the Hamiltonian 
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where Ci (cj) is the annihilation (creation) operator 
at site i, rii = c\ci the number operator, and the 
flux-dependent boundary condition enters through cq = 
exp(z<I>)ci. The interaction is restricted to nearest neigh- 
bors and effective in the sample, but vanishing in the 
lead. It is equilibrated by a compensating potential 
that prevents the particles from emptying the interact- 
ing region. The form of the Hamiltonian allows to have 
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FIG. 2: Conductance g as a function of the interaction 
strength U for different values of the sample length Ls- 

particle-hole symmetry at half filling. We work with a 
number of fermions N = L/2, such that the mean den- 
sity is always 1/2 independently of Lg and Ll. 

Using the DMRG algorithm as described in Ref. [ p^ , 
we calculate the ground-state energies £'($) at $ = 
and $ = TT, to obtain the stiffness D — {L/2)\E{Q) — 
i?(7r)| (which is simpler to calculate than J). When the 
flux-dependence of the ground-state energy is sinusoidal, 
which is the case in one-dimensional localized systems, D 
is directly proportional to the maximum J. We will use 
such an identification to calculate g from Eq. , taking 
D instead of J |||. 

An obvious requirement for Eq. to be useful is that 
the resulting g should be independent of Ll, since the 
lead is an auxiliary element of the procedure. Therefore, 
the first numerical step is to compute D for increasing 
Xl with given ig and U (Fig. We find a very clear 
asymptotic behavior, described by the scaling law 

D{U,Ls,U)^D^{U,Ls) e^^{^^^y (7) 

where the intrinsic value Daa is independent of the length 
of the auxiliary lead. This asymptotic value Dao is then 
used to determine the conductance as g = {D^/ , 
where corresponds to the clean non-interacting ring. 
The sign of C depends on the parity of the number of 
sites Lg in the sample: C > for even Lg (filled symbols) 
and C < for odd Lg (open symbols). 

The method also works for non-polarized electrons 
(Hubbard model with on-site interaction), and the size- 
scaling again allows to obtain intrinsic values (lower inset 
in Fig. |l|). In the sequel we concentrate on the spinless 
case (Hamiltonian (^)), which contains the main features 
we are interested in, and allows to reach larger samples. 

Having verified the consistency of our approach, we 
now study the systems of interest. In Fig. we present 



FIG. 3: Conductance g as a function of the interaction 
strength U , for a fixed L^s — 12, and increasing smoothing 
of the contacts, defined by the length La (see inset). Filled 
(open) symbols correspond to even (odd) Ls; La ~ 1 (o), 
and La = 9 (□), La = (■), La = 4 (•), and La = 10 (a). 
Using the same smoothing length (La = 10) but improving 
in the shape (a tanh function (thick solid line) instead of a 
linear increase (a)) helps g to approach the perfect value. 



the conductance, as a function of U, for various sam- 
ple lengths Lg. One observes a monotonic decrease with 
increasing U , and a very clear even-odd asymmetry ac- 
cording to the parity of Lg. Samples with odd Lg ex- 
hibit almost perfect transmission up to large values of U, 
while an even Lg results in a decrease of g{U) already 
at weak interaction. For odd Lg , particle hole symmetry 
leads to degenerate sample configurations with (Lg±l)/2 
particles in the interacting region. This is similar to a 
Coulomb blockade resonance. The traveling particle can 
thus become trapped for a long time (tw > 0), consis- 
tent with negative C Since the two configurations 
are coupled by processes which transfer particles through 
the interacting sample, one obtains perfect transmission, 
as in the perturbative treatment of Ref. . 

On the other hand, an even number of sites implies 
that the transport of one particle through the sample 
takes place via a virtual state with an energy of order U 
above the ground-state. Thus, no resonance can be ex- 
pected and the transmission, which is suppressed already 
by moderate f7, is a fast process with tw < 0, consistent 
with C > 0. In addition, increasing Lg reduces g linearly 
for small U, and exponentially for U > 2, consistent with 
the Mott-insulating behavior. 

The even-odd asymmetry, and the perfect transmission 
for the odd case, point to the importance of the contacts. 
In order to investigate their role, we introduce a posi- 
tion dependent interaction strength Ui which increases 
linearly from to its maximum value U, inside the "con- 
tacts" of length La (see inset of Fig. As shown in 



4 








1 2 3 4 5 6 

U 

FIG. 4: Ensemble averages of In <; as a function of the inter- 
action at Ls ~ 8, in the presence of disorder W = 1 (a), 
5 (•) and 9 (♦). The statistical errors are smaller than the 
symbol size. The dash-dotted line represents the clean case, 
the dotted Unes represent four individual samples ai W = 9. 



moderate repulsive U, showing the non-trivial interplay 
of disorder and interactions in a transport problem. 

In conclusion, starting from a recent proposal JlO| , we 
have provided a well-defined procedure for calculating 
the conductance g of interacting one-dimensional wires, 
and used it to investigate correlation and disorder effects. 
While the interaction reduces g for spinless fermions in 
the presence of weak or moderate disorder, a moderate re- 
pulsive interaction increases g at strong disorder. We also 
determined the crucial role of the sample-to-lead contacts 
on the conductance. 
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programs. RAM acknowledges the financial support from 
the European Union's Human Potential Program (con- 
tract HPRN-CT-2000-00144). RAJ and DW thank the 
INT at the University of Washington for its hospitality 
and support during completion of this work. 



Fig. |[ these smooth contacts increase g, at constant ef- 
fective length Lcs — — Lp^ of the sample. The effect 
is pronounced in the case of even Lg and the conduc- 
tance approaches the ideal situation of perfect transmis- 
sion expected for rcflcctionless contacts when we improve 
the smoothing While the form of the contacts pre- 
serves the right-left symmetry, the smoothing leads to an 
extension of the perfect transmission at odd Ls towards 
higher values of U . However, we found that asymmet- 
ric contacts destroy the perfect transmission. The strong 
influence of the contacts is crucial when describing exper- 
iments since it seems impossible to connect a nanosample 
via reflectionless contacts. It also shows the limitation of 
other approaches relating the conductance of an interact- 
ing sample to its intrinsic properties, without taking into 
account the way it is connected to the electrodes. 

While clean interacting systems are of physical interest 
{e.g. carbon nanotubes), it is also important to consider 
the generic case of disordered systems. To this end, we 
add the term W X^i^i ''^i'^i to the Hamiltonian (^), where 
W is the disorder strength and the vi are distributed 
equally in [—1/2,1/2]. We have checked that D scales 
with L as before (Eq. (Q)), ensuring a well-defined limit- 
ing value for g. The even-odd dependence of g disappears 
when disorder is introduced. The combined effect of dis- 
order and interactions on g is shown in Fig. ^ In the 
ballistic regime (at = 1 the mean free path exceeds 
Ls) the effect of W is weak at small J7, and it becomes 
more pronounced at stronger U (when the disorder pins 
the Mott insulator, reinforcing the localization). At large 
W , g for individual samples exhibits peaks as a function 
of U whenever a charge reorganization occurs, similarly 
to the case of the D calculated without the lead |2^ . Very 
remarkably, the ensemble average of In g is increased by a 
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